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An exact solution in the Faxen approximation is given for the Lamm equation in which the sedi- 
mentation coefficient is related to concentration as s=$>(l — kc). It is shown that the solution in this 
case can be expressed in terms of the solution to the linear case (k=0) with a modified argument. 
The boundary sharpening phenomenon expresses itself very clearly in the solution presented here. 
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Solutions of the Lamm equation for constant diffusion and sedimentation coefficients have 
found wide application in the analysis of data from the ultracentrifuge [l]. 1 It is known that both 
the sedimentation and the diffusion coefficients depend on concentration, and in particular, the 
sedimentation coefficient can be approximated by 



s = s l(l + kc) (1) 

where s is the sedimentation coefficient at infinite dilution and k is a constant which is a measure 
of nonlinearity. Fujita has shown that if the relation of eq (1) is replaced by 

s = s (l-kc) (2) 

while the diffusion coefficient remains independent of concentration, the resulting Lamm equation 
can be linearized rigorously [2]. Fujita solved this Lamm equation in the Faxen regime by an 
approximation technique. It is the purpose of this paper to show that an exact solution of the 
Lamm equation is possible in the Faxen regime, and that the solution for the nonlinear case can be 
written entirely in terms of the solution to the linear problem. The technique exploited in this 
paper was originally suggested in a preliminary note [3] although the result given there contains 
an error. Our solution then enables one to derive an expansion in powers of the time variable 
to as high a degree as desired, whereas Fujita's original solution was limited in this respect. The 
Lamm equation can also be solved with the ansatz of eq (2) in the rectangular approximation [4], 
and an Archibald solution which takes account of boundary conditions can be written [5]. A de- 
tailed comparison between solutions of the Lamm equation with the two forms of sedimentation 
coefficient given by eqs (1) and (2) will be presented elsewhere [6]. 
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The Lamm equation will be written 

dc ldZ-.dc „ 0/ ., , v \ 

— = -- rD — -sotoVl-Acc (3) 

dt rdr\ dr ] ; 

in terms of the solute concentration c, the diffusion coefficient D, sedimentation coefficient so, and 
angular frequency a). Let r be the radial position of the meniscus, and Co be the initial concen- 
tration. Then the transformations 

0=clco, e = 2D/(so(x) 2 rf) x = (r/r ) 2 , T = 2a) 2 s t a = kc (4) 

convert eq (3) into 

This equation can be linearized by introducing a new dependent variable u(x, r) by 

0(*, t)=- — In u(x, r). (6) 

It proves convenient to also define new independent variables by 

z = 2(xe- T y' 2 , £=l-e- T (7) 

which, together with eqs (5) and (6) implies that u satisfies 



1 du _ d 2 a 1 du 

€ d£ dz 2 z dz 



(8) 



which is to be solved with the initial condition 

u{(T, 0)=1, 0^o-^2 = exp^(cr 2 -4), a > 2. (9) 

We can solve eq (8) by a separation of variables. A general solution for which dujdz remains 
finite at x = is 



a(z, £) = ^ + £z 2 + z | X g(X)/i(Xz)c- x ^rfX 

Jo 



(10) 



where ^ and 5 are constants, and g{\) is to be determined from initial conditions. Notice that the 
first two terms are solutions to the homogeneous equation obtained by setting dzz/d£ = 0. Initially 
we have 

u(z, 0)-A-Bz 2 = z f°° g(k)J l (Xz)dk (11) 

where g(X) can be calculated by regarding the right-hand side as a Hankel transform [6]. In this 
way we find 

g(X) = X | o °° 7 i(Xcr)Nc7, 0)-A- Ba 2 ]da. (12) 
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Hence u(z, £) is given by [7] 

u(z l = A + Bz 2 + z I™ [u(o-,0)-A-B(T 2 \d(T f °° AJi(XcrX/i(Az)e- x ^d\ 
Jo J o 



= A+Bz 2 + — e-^ 



/o" C ~^ A (§) M °"' )"^~ Bcr2 ] dcr < 13 ) 



where Ii(/jl) is a Bessel function of the first kind, with imaginary argument. 
Since 



J o- 2 e-^ A ^j do- = 2ze£e^ 



(14) 



we see that the terms containing B cancel. The constant A can be determined by requiring that 
the limiting solution for a — » agree with the known Faxen solution for constant s and D. It can 
easily be verified that this is equivalent to setting A = l. Hence an expression for u(z, £) is 



*. 0- 1+gJ/ «"*=?/. (g) [«*«-« L l]*r. 



(15) 



It will prove convenient for later purposes to perform an integration by parts, noting that I\(y) 
— dIo{y)/dy. This leads to the expression 



u(z, 0= 1 +^/ 2 ' o-e-^o (g) [(l-o0e5<"*-«- l]«ftr. 



(16) 



Differentiation according to eq (6) together with an integration by parts leads to an expression* 
for 0(z, 0, 






(17) 



When a= this relation reduces to the known Faxen solution, which will be denoted by #o(z, £): 

/w y\ i — £ f x (<r 2 + z 2 ) T I Z(T\ , 

We now notice that all of the integrals appearing in the expression for 0(z, £) have the same general 
form of an exponential of a quadratic in a, multiplied by the common factor 



O" exp 



(a 2 + z 2 ) 



4e£ 



h 



zcr 
2e£ 



Hence all of the integrals are expressible in terms of O provided that the arguments are properly 
chosen. After some algebra it is found that 



0(z, o 



_ l-{ 



e—y-n^Mz*,? 



l-(q+l){ l-q£ « i f_ 1 „ „ 

1 + 1 / _i_ i\y e e 1 4ll-aC)JCo(Z , t 

1 — (a + Ik 



(19) 



1-C 



6>o(z, 
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where 

i^rS'i^i^T (20) 

1 — a£ 1 — at, v ' 

Simple properties of the solution follow immediately from the representation of eq (19). The 
parameter a/e which appears in the expression of eq (19) is commonly greater than 10 for values of 
D, 5, and o) 2 (which validate the approximation of eq (2)) and the Faxen approximation. Hence 
the exponential terms in eq (19) predominate when 

1 (21) 



4(1 -oQ 
is positive by an order of magnitude greater than e/a. When this is the case we have 

6(z, Q~-^^ = [a + (l- a)er i (22) 

which defines the behavior of the plateau region. This can also be verified directly from eq (5). 
When the quantity of eq (21) is negative by an order of magnitude greater than e/a, d(z, £) is effec- 
tively equal to zero. The transition region occurs from values of r for which 

1=0 (-,)■ 



4(1 -aQ 



Thus, in the limit e = the transition region (from essentially zero to a plateau value) is very sharp, 
and for fixed e and increasing a the transition region also becomes sharper, as indeed is well known. 
If the gradients are denoted by G(z, Q and Go(z, £), i.e., 

then we can write for G(z, £): 

c , 1 G4?,Q | l-^{ l -lh MzA \ 

[,u l-a£F(z,z*, C, {*) l-a{ F 2 (z, z*, £, £*) 

x ter^TTK-^'-^^*' ^ + Tifew G ^*' D-^GoU, o} (25) 

where 

F(z,z*, C, £*)=! + 1 }~f lK e-- € [ l -W^^ (z*, H-JZ^Wz, 0- (26) 

It is possible now to derive successive approximations to 6(z, Q by expanding it in powers of 
(e£) 1/2 . That this is the natural parameter for the expansion will be seen later. Since 6(z, Q 
depends only on do(z, Q we must first consider the expansion of the latter function for small £. 

For small £ the Bessel function /o \w~p) nas the asymptotic expansion [7] 
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/here 



\2e£/ Vttzo- j^ 



(-iy IX/+D 



Therefore 6*0(2, £) has the expansion 
e~ 



0<>(z, - 



2v*ze£jtj 



2Bj ,MYf%-V^ 



7-1/2 



fc|M¥)7; 



(27) 



(28) 



-l|2/2 



dw. 



(z+i/.v2€jy- i/2 



(29) 



We now consider the integrals which appear in this series and expand them in powers of (e£) 1/2 . 
To carry out this expansion we write, for fixed u and small £ 



r=o \ r / 
This can be further transformed by noting the identity 

(-J+ll2\ ( _ 1Y (r + j-3l2\ 



(30) 



(31) 



Finally, combining eqs (29), (30), and (31) we find that the complete asymptotic expansion of the 
Faxen solution can be written as 



»*.0~SitHr S ,(»-'-f)^»„fi 

V 27T j = ,n = 2j \ ./-3/2 / Z'" W, 

where the U r (y) are defined by 

UAy)= J * We-* 1 * 



(III. 



In particular, the first two terms of eq (32) yield the usual Faxen approximation 

2-z\ , 1 (et\ 11 * 



V^r I \V2€£/ z \2> 



D«v*] 



(32) 



(33) 



(34) 



in which Udy) can be expressed in terms of an error function. A complete asymptotic expansion 
for 0(z, Q can now De obtained by combining eqs (19) and (32). Similarly, 30/d(r/ro) is expressible 
in terms of 0o(z, Q and its derivative, Go(z, J), which has the expansion 



V27T j=0 m=2j \ J 6/Z / Z 



• Gg) 



-fc/,„_, f^U- (2 



2 C {/ (2e£)2t§£! e ' * { 



(35) 



The lowest order Faxen approximation to 0(z, £) is 



0(z, £) = 0„(£) 



!-</> 



1- 



2-z \ 
V2e£/ 



^^•['-^H^ 



(36) 
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where P (Q is the normalized plateau concentration (eq (22)) and (/>(/*) = (27r)~ 1/2 exp (—x 2 /2) 

Jix 

dx is the complementary error function. To this approximation, the gradient is 



36 n l9 d0 e- Tl2 ~^[ l ~W^o\ 



d(rlro) dz D 



iTT^f^K^+TT^fl-^MlW 



rl / I ^ 4 ' t b ' / \ -*- rl / 

7re£ \V2eT/ V^£* \ \V2^ 



e(l-aQ \V2e£*/ V \V2e£ 



:*lfs=)li-*(r^=)H (37) 



fhere 



2-fV 



When a is set equal to eq (37) reduces to the known Faxen solution. 

Fujita [2] showed that in his approximation the position of the maximum is given by 

r max =r e T [l-oi(l-e-% (39) 

The result of eq (19) implies that this is modified to 

c =r o e'[l-a(l-e-0][l+/]. (40) 



f'n 



The factor /is generally of the order of 0.1 or less and therefore Fujita's result for the position of 
the maximum gradient is sufficient for all practical purposes. It also follows that Fujita's calcula- 
tions for the height-area ratio are sufficiently good for practical purposes. 
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